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We consider theoretically the effects of an applied uniform magnetic field on the magnetic spec- 
trum of anisotropic two-dimensional and Dzyaloshinskii-Moriya layered quantum Heisenberg anti- 
ferromagnets. The former case is relevant for systems such as the two-dimensional square lattice 
antiferromagnet S^CuC^Cb, while the latter is known to be relevant to the physics of the layered 
orthorhombic antiferromagnet La2Cu04. We first establish the correspondence between the low- 
energy spectrum obtained within the anisotropic non-linear sigma model and by means of the spin- 
wave approximation for a standard easy-axis antiferromagent. Then, we focus on the field-theory 
approach to calculate the magnetic-field dependence of the magnon gaps and spectral intensities for 
magnetic fields applied along the three possible crystallographic directions. We discuss the various 
possible ground states and their evolution with temperature for the different field orientations, and 
the occurrence of spin-flop transitions for fields perpendicular to the layers (transverse fields) as 
well as for fields along the easy axis (longitudinal fields). Measurements of the one-magnon Raman 
spectrum in S^CuChCh and La2Cu04 and a comparison between the experimental results and the 
predictions of the present theory will be reported in part II of this research work [L. Benfatto et al, 
forthcoming article, cond-mat/0602664 . 
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I. INTRODUCTION 

The relevance of antisymmetric superexchange interactions in spin Hamiltonians describing quantum antiferromag- 
netic (AF) systems has been acknowledged long ago by DzyaloshinskiiiA Soon after, Moriya showed that such interac- 
tions arise naturally in perturbation theory due to the spin-orbit coupling in magnetic systems with low symmetry? 2 . 
Nowadays, a number of AF systems are known to belong to the class of Dzyaloshinskii-Moriya (DM) magnets, such 
as a-Fe203, LaMn03y2> and P^VaOsr 1 and to exhibit unusual and interesting magnetic properties in the presence of 
quantum fluctuations and/or applied magnetic fieldi ^'^ 

Also belonging to the class of DM antiferromagnets is La2Cu(J4, which is a parent compound of high-temperature 
superconductors. In La2CuC>4 the unique combination of antisymmetric superexchange, caused by the staggered tilting 
pattern of oxygen octahedra around each copper ion in the low-temperature orthorhombic (LTO) phase, and weak 
interlayer coupling, results in an interesting four sublattice structure for the antiferromagnetism of La2Cu(J4, where 
the Cu ++ spins are canted out of the CuC>2 layers with opposite canting directions between neighboring layers JiMi 1 ^ 
These small ferromagnetic moments lead to a quite unconventional physics in the antiferromagentic phase. For 
example, when a small magnetic field is applied perpendicular to the layers the magnetic susceptibility shows a strong 
enhancement as the Neel temperature is approached, due to the formation of these ferromagnetic momentsi ^ 12 ' 1 ^ . 
When the field strenght is further enhanced one can observe a spin flop of the ferromagnetic moments whith respect 
to the out-of-plane staggered order that they display at zero field&Afiiii Analogously the spin-flop transition of the 
staggered in-plane AF moments for a field along the easy axis is accompanied by new effects related to the presence of 
the DM interaction. 9 " 10 Recently new theoretical approach es hav e been propose d 12 ' 1 ^ 14 to integrate the semiclassical 
picture already presented in the prior work of Refs. [?Hftl lflllOl |. Even thought the basic physical picture remains 
the sa me, the inclusion of quantum effects in the long-wavelenght formulation of the spin problem discussed in Ref. 

allowed for a straightforward and complete understanding of the unusual magnetic-susceptibility anisotropics 
observed in La^CuO^ for a rather large temperature range, < T < 400 K. Moreover, a particular attention has 
been devoted in Ref. [lj| to the analysis of the one-magnon excitations by means of Raman spectroscopy, and the use 
of the long-wavelength theory turns out to be very convenient to understand why the DM interaction is behind the 
appearance 15 of a field-induced mode for an in-plane magnetic fields 

A better understanding of the anomalies related to the presence of the DM interaction in La2Cu(J4 compounds 
can be achieved by directly comparing its properties with those of a similar spin system like S^CuCbC^. In this 
case the DM interaction is absent due to the higher crystal symmetry, but spin-orbit coupling can still give rise to 
small anisotropics of purely quantum mechanical origin^ Thus quantum corrections coming from spin-orbit coupling 
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give rise to a quite small easy-axis anisotropy, so that a gap in the in-plane magnon excitations has been observed 
in ESR. 17 As a consequence, S^CuCbC^ behaves as an ordinary easy-axis antiferromagnet, in contrast to La2CuC>4 
which should be classified as an unconventional easy- axis antiferromagnet. 

It is the purpose of this article to study in detail the influence of an applied uniform magnetic field on the magnetic 
spectrum for the above two cases: anisotropic two-dimensional and layered DM antiferromagnets, and to compare 
the qualitative differences between these two cases. This will be done mainly using the continuum quantum field 
theory appropriate for each of these two cases, i.e. the non-linear sigma model (NLSM), properly modified to account 
for conventional or DM anisotropies. Nonetheless, we will sketch in the beginning the calculation of the magnon 
gaps within the framework of the semiclassical approximation for conventional antiferromagnets, and we demonstrate 
the complete equivalence between the two approaches as far as the gaps values at low temperature are concerned. 
However, as it will become clear in the following, the quantum NLSM followed here allows also to account for the 
quantum and thermal effects of the spin fluctuation, which were neglected in the previous approaches&Siifl. As a 
consequence, we can evaluate (within the given saddle-point approximation for the transverse spin fluctuations) the 
full (H, T) phase diagram for a field in the various direction, which can be compared with the existing experimental 
data. At the same time, the continuum field theory provides an elegant and straighforward description of the spin 
fluctuations in the coupled-layers case, and also of the various spin-flop transitions that may occur in La2CuC>4 at 
moderate fields. 

The structure of the paper is as follows. In Section II we introduce the model Hamiltonians for Sr 2 Cu02Cl2 and 
La2CuC>4, appropriate to describe a conventional anisotropic two-dimensional antiferromagnet and a DM antiferro- 
magnet, respectively. In Section III we sketch the standard semiclassical calculations of the magnon gaps at zero 
field in the two cases, and at finite magnetic field for the standard anisotropic case. Then the same (conventional) 
results are reproduced in Section IV using the NLSM approach, whose main properties are here described. Section V 
is dedicated to the case of a layered DM antiferromagnet, and the effects of an uniform magnetic field applied along 
the three crystallographic directions are extensively discussed, with reference to the specific structure of La2CuC>4. 
The conclusions are reported in Section VI. In a second partji& we shall make a quantitative comparison between 
the predictions of the theory developed in this article and the magnetic spectrum probed by one-magnon Raman 
scattering in both La2Cu04 and Sr2Cu02Cl2- 



II. SPIN-HAMILTONIANS FOR CONVENTIONAL AND DM ANISOTROPIES 



La2CuC>4 is a body centered orthorhombic antiferromagnet with Bmab crystal structure. A single layer of CuC>2 
ions in La2CuC>4 can be described by the S = 1/2 Hamiltonian 

H sl [S, D] = J ^ • + ]T D 4j • (Si x Si) + s i • r, 3 • Sj , (1) 

where and T y are, respectively, the DM and XY anisotropic interaction terms that arise due to the spin-orbit 
coupling and direct-exchange in the low-temperature orthorhombic (LTO) phase of La2Cu04i£ Throughout this work 
we adopt the LTO (abc) coordinate system of Fig. for both the spin and lattice degrees of freedom, and we use 
units where % = kg = 1. 

The direction and the alternating pattern of the DM vectors, shown in Fig. ^ have been calculated by several 
authors^ by taking into account the tilting structure of the oxygen octahedra and of the symmetry of the La2CuC>4 
crystal. For La2Cu04 the DM vectors are in good approximation perpendicular to the Cu — Cu bonds and change 
sign from one bond to the next one: 

B AB = ^=(d,-d,0), B AC = ^=(d,d,0), (2) 
while the XY matrices T provide essentially an easy-plane anisotropy for the Hamiltonian Q: 

r ab= 

where AB and AC label the Cu ++ sites on horizontal/ vertical bonds respectively (see Fig. As it has been 
stressed by Shekhtman et alii, even though the parameters d and ^,2,3 > have different orders of magnitude, 
with d r~j 10~ 2 J and Ti ~ 10 _4 J, they should be considered on the same footing (see also discussion following Eq. 
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Ill2t below). Indeed, one can show that considering the two last terms of Eq. the interaction between spins on a 
neighboring bond can be written in a completely isotropic form by rotating locally the spin operators around the Dy 
axis by an angle % = arctan |D,j |/2 J. As a consequence, weak ferromagnetism arises only when global frustration 
of the DM pattern exists. In term of the DM vectors defined above, this condition corresponds to having d + ^ d_, 
where d± = (Dab ± Dbc)/2. This condition is clearly satisfied by the DM vectors Q, so that WF is expected in 
La2Cu04- 

A realistic model for La2CuC>4 should include also interlayer coupling. In the orthorhombic unit cell of La2Cu04 
the spins of the Cu atoms are displaced by an in-plane diagonal vector (1/2, 1/2, 0) from one layer to the next one. 
As a consequence, given a couple of spins in a layer and the nearest couple in the next one, the DM vector of the 
corresponding bond will change sign. We can then write the full Hamiltonian as: 

H = J ± S m ' S m+1 + H *l I 13 ™ - ° m ] . (3) 

m m 

where S m represents the spin at a generic position of the rath plane and D™ B Ac = (—l) m T)AB,AC ■ Since the 
La2CuC>4 unit cell is body centered, the coupling J± in Eq. (j^J connects the two spins at (0,0,0) and (1/2,0,1/2) 
in the LTO notation. It is worth noting that the pure 2D system does not display any spin rotational symmetry, 
so it can order at finite temperature without violating the Mermin- Wagner theorem. However, in the presence of an 
interlayer coupling the transition to the AF state will ultimately have 3D character, with spins aligned AF also in the 
c direction. On the other hand, as we shall discuss below, the interplay between the interlayer coupling and the DM 
interaction can lead to a quite unconventional behavior in the presence of a finite magnetic field. Indeed, the DM 
interaction is not only the source of an easy-plane anisotropy (the spins prefers to align perpendicularly to the DM 
vector d+), but induces also an anomalous coupling between the AF order parameter and an applied magnetic field. 

These effects can be better understood by comparing the results obtained with the model (and its three- 
dimensional version @) with the ones coming from a more conventional 2D anisotropic Heisenberg model, as: 

H con = JS i S *j + ( J - a a )S?S^ + ( J - a c )5, c 5j. (4) 

(i,i) 

The Hamiltonian l@J is the appropriate starting model for SraCuOaCla, where interlayer coupling is even less relevant 
than in La2CuC>4 due to the frustation on the tetragonal unit cell. Here the crystallographic in-plane a, b axes are 
choosen with b parallel to the spin easy axis (which is along the xy direction) , so that we will have a similar notation 
to the one used for La2Cu04. However, a = b for Sr2Cu02Cl2, since the system is tetragonal. As we explained in 
the introduction, a a should be zero in a tetragonal system. Nonetheless, quantum effects can induce an in-plane 
anisotropy 1 ^ which we will mimic with a finite a a anisotropy term in what follows. To clarify to what extent the DM 
interaction introduces an anomalous behavior, we shall start our analysis of easy-axis antiferromagnetism from the 
anisotropic model <@}. We will then be able to go back to the model an d to correctly distinguish the effects of 

the magnetic field alone from the ones arising from the presence of the DM interaction. 




FIG. 1: Left: the hatched circles represent the O ions tilted above the CuCh plane; the empty ones are tilted below it; small 
black circles are Cu ++ ions. Right: Schematic arrangement of the staggered magnetization (small arrows) and DM vectors 
(open arrows). Right bottom: definition of the vector d+ = D+/4S 1 . 
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III. SEMICLASSICAL APPROACH 



A. Conventional anisotropics 

Let us start our analysis of the single-layer Hamiltonian using a semiclassical approach. Similar calculations 
have been already carried out in different contexts j£L2L and here we shall just review the main steps to fix the notation 
and to clarify the conventional behavior of an ordinary easy-axis antiferromagnet. Let us denote by Si and S2 the 
spins on the two AF sublattices. The free energy density in the AF phase can be written: 

F = z JSfS$ + z(J - a a )S?S2 + z(J - a c )SlS c 2l (5) 

where z is the number of nearest neighbors and S a ' b,c are the components of the vector along the three crystallographic 
axes. Within the semiclassical approach, the spins are treated as classical vectors of length S: thus the ground-state 
configuration can be easily determined by imposing that Vs^-F = 0. This condition clearly shows that the spins order 
along the b direction, with S° = — Sf] = Sx b . To calculate the magnon gaps one uses the classical equations of motion: 



dSj 
~dt 



S, x V Sl F, 



1,2 



(6) 



where Vs^/ls^Ms) represents the effective local magnetic field around which each magnetic moment (<7 s //b)S pre- 
cesses. Here g s is the gyromagnetic ratio and \ib is the Bohr magneton. By expanding the spins around the classical 
solution, Si = S(m®, l,mj), S2 = S^mf,— 1,"4) Eqs. © give a set of 4 coupled equations for the time dependence 
of the transverse spin fluctuations, which can be easily solved by putting mf(t) — mf exp(iuit). It then follows that 
the ui are the eigenvalues of the matrix (corresponding to the vector (mf , m|, mf , mf)): 
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(7) 



The eigenvalues uj = ±uj a and u = ±uj c correspond to eigenmodes describing spin fluctuations with a larger a/c 
component respectively, allowing for the identification of uj a and uj c as the magnon gaps for the in-plane and out-of- 
plane spin- wave modest 



zSy/2a a J, 



zS\/2ol c J \ 



(8) 



where we approximated y / 2a a , c J oi a a c « y / 2a a , c J because a ai a c <C J. 

In the presence of a finite magnetic field a term — J2i H • Sj must be added to Eq. J3J, which translates into a term 
— H ■ (Si + S2) in Eq. ©. Observe that in what follows we shall measure the magnetic field in units of g s ^B = 1, 
unless explicitly stated. Then one follows the same procedure as before, by noticing that when a transverse field is 
applied (i.e. a field perpendicular to the easy axis) the sublattice ground-state configurations S? 2 acquire an uniform 
component in the field direction proportional to H/z(2J — a a , c ) (for a field along a and c, respectively). By adding 
fluctuations transverse to the new equilibrium direction one finds for example for H || a the fluctuation matrix: 



zS 



( J J-a c \ 

-J + a c -J 

-J cos 2 + (J - a) sin 2 (j) - (H/zS) sin 4> J sin 2 (f> — (J — a a ) cos 2 <j) 

y - J sin 2 <t> + (J - a) cos 2 J cos 2 </> - (J - a) sin 2 <p + (H/zS) sin <p J 



where 4> = arcsin(i? / zS(2J — a a )) ~ H/zS2J. As a consequence, the new magnon gaps areiSl 

u) a (H) = y'ml + H 2 , uj c (H) = m c . 



(9) 



Observe that since the canting of the spins due to the magnetic field is small, the eigenmodes still describe fluctuations 
having predominantly a or c character, respectively. We see that the effect of a transverse magnetic field is to harden 
the gap of the mode in the field direction, and to leave the other gap unchanged. Indeed, when H is parallel to c we 
find a similar result, with an increasing ui c gap and a constant ui a gap. 



5 



Finally, let us consider the case of a longitudinal field, i.e. of a field parallel to the easy axis. In this configuration 
no uniform spin magnetization develops, but the magnetic field effectively shifts the AF coupling along the easy axis 
in the two sublattices, so that the new fluctuation matrix reads: 



zS 



( J + H/zS J-a c 

-J + a c -(J-H/zS) 

-(J + H/zS) -J + a a 

V J-a a J-H/zS 



The four eigenvalues are given by 



lu 2 /(zS) 2 = J(a a + a c ) - a a a c + {H/zSf ± ^{a c - a a ) 2 J 2 + 4(a a + a c )J(H/zS) 2 - (H/zS) 2 {a a + a c ) 2 , 
and using the fact that a a , a c <§; J they can be readily expressed in terms of the bare gaps m a , m c as: 



m 2 + m 2 



H 2 



+ H 2 



AH 2 



m 2 + m 2 



AH 2 



(10) 



where we assumed m c > m a . Note that at fields larger than the bare gaps one observes essentially a linear increase of 
the magnon gaps with the magnetic field. In the case of degenerate gaps, a c = a a , only the linear regime is accessible 
and Eqs. JEll simplify to&L 



H, u c 



H. 



B. Dzyaloshinskii-Moriya interactions 

Let us discuss how the previous results are modified in the presence of DM interactions. First, taking into account 
that ri j2 <C d <§; J the free energy density corresponding to the Hamiltonian can be written as: 

F = z.J{SlS% + S^) + z(J - a c )SlS L 2 + zd+{S\S c 2 - SfS%), 

where a c = Ti — > and d + = |d+| = d/y2. The ground-state configuration of the previous free energy has been 
already discussed by several authors (see for exampl o 7 i 8 i 10 ). The spins order AF with an easy axis in the b direction, but 
with an additional small ferromagnetic (FM) component along c, = 5(0, cos cj>, sin </>) and S° = 5(0, — cos<^>, sine/)). 
The angle (j> of the out-of-plane canting of the spins is given by 0o = (1/2) arctan[2rf+/(2J — at c )] ~ d + /2J and it is 
due to the DM interaction (see also Fig. below). When this canting is taken into account in the linearized equations 
of motion © , one can easily see that the matrix for the transverse fluctuations in zero field has the same structure 
of Eq. Q, with 

a a = d+<l>o w d\/{2J). (11) 

As a consequence, the effect of the DM interaction is twofold: it induces the FM canting of the spins, and it reduces 
the AF coupling in the a direction. The corresponding magnon gaps are, using Eq. (|HJ and the equivalence (tTTll : 

uj a = m a = zSd+, lu c = m c = zS\/la c 3 . (12) 

Notice that the gap of the a mode is proportional to d, while the gap of the c mode scales with the square-root 
of the parameter a c = T% — T3. As a consequence, even though Ti ~ 10~ 4 J and d ~ 10~ 2 J the two gaps are of 
the same order of magnitude in La2CuC>4. When a finite magnetic field is applied the system will evolve towards a 
new ground-state configuration. Following the procedure described above the new magnon gaps can be determined. 
However, we will not present these calculations here, because we shall describe in detail in the next sections how these 
results can be obtained using the NLSM approach, and how do they differ from the results Q and 11011 . that we shall 
refer to as 'conventional' in what follows. 
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IV. NLSM FORMULATION FOR CONVENTIONAL ANISOTROPIES 

In this section we will show how the behavior of the spin gaps in the presence of magnetic field can be easily derived 
within a NLSM description of the low-energy physics of the spin model (@J and ©■ We will first discuss the simple 
anisotropic model (QJ, to show the agreement with the results © and iflOj) presented above. Since the semiclassical 
approach is much more lengthly and less transparent than the NLSM description, we shall adopt the latter to deal 
with the more complicated case of the Hamiltonian Q. 

The derivation of the NLSM starting for the 2D Heisenberg Hamiltonian has been extensively discussed in the 
literature^ Here we just recall the main steps and stress the origin of the mass terms due to the anisotropies in 
Eq. First, we decompose the unit vector fli = S;/5 at site Tj into its slowly- varying staggered and uniform 
components, 

^ = = eiQ ' X,n ( x - r ) + aL ( x - r )' ( 13 ) 

where Q = (ir, n) and a is the lattice parameter. The constraint f2f = 1 is enforced by nf — 1 and Li • = 0. Using 
this decomposition, the Heisenberg part of the Hamiltonian Q has the standard formal 

TC-2 r 

Chj = J J2 S * ' S i = — / d2x t( Vn ) 2 + 8L2 ] . ( 14 ) 
while the terms proportional to a a , a c give rise toi 2 ^ 



S 2 f 

C a = — aaS^Sj + a c S^S C j = — / d 2 j 



2a a n 2 a + 2a c n 2 c + -^a 2 a a (\7n Q ) 2 - ^a 2 a c (\7n c f 



(15) 



Since a a<c <C J we can neglect the correction induced by the small anisotropies to the gradient of the transverse modes 
in Eq. l|Ulj) . and we can retain just the first two terms of L a . Using a path- integral coherent states representation of 
the spin states, which in addition to the previous contributions gives rise to the (dynamical) Wess-Zumino term 2 ^ 



— iS f 

Cwz = / d 2 xL • (n x n), 

a J 



we can obtain the partition function Z = J DnS(n 2 — l)e 5 , with the action S = J At[Lhj + £ a + Cwz]- After inte- 
gration of the L fluctuations we obtain the following anisotropic non- linear a model (/3 = 1/T and J — J Q dr J d 2 x): 



So = ^ J {(drn) 2 + c 2 (Vn) 2 + m 2 n 2 + m 2 n 2 } . (16) 

The bare coupling constant g and spin velocity c are given by gc = 8 Ja 2 and c = 2y/2JSa, and we defined m 2 c — 
32JS 2 a ac , which corresponds to the result |JH} above with z = 4, as appropriate in two dimensions. In generic d 
dimensions the coefficients 2a a , c in Eq. ltl5|l are replaced by z/2a a , c and one puts gc = 4:dJa d — 2zJa d , leading again 
to the definition ® of the masses. In the NLSM itlfill the spin stiffness is renormalized by quantum fluctuations to 
p s = c(l/A r g — A/47r)(22i 2 ^i2& where A is a cutoff for momentum integrals and N — 3 is the number of spin components. 
When the system orders we find the staggered magnetization at (n) = onx;,, because the orientation along x a or x c 
would cost an energy m a or m c , respectively. 

It is worth noting that even though the NLSM llKit contains explicitly only the staggered spin component n, 
nonetheless the saddle-point value of the uniform spin component L determined before integrating it out contains the 
residual information about the ferromagnetically ordered spin component. This is evident when an external magnetic 
field is applied on the system. In this case one can easily repeat the previous calculations by taking into account that 
the saddle-point value of the uniform magnetization L acquires an additional contribution proportional to Hi2£ 

; (nxn) + -i-[H-n(n-H)]. (17) 



8aSJ" ' 8aSJ l 

Observe that the first term is proportional to the time derivative of n, so it averages to zero for the equlibrium 
configuration (indeed no FM component is present in the ordinary AF phase). However, at finite H a non- vanishing 
average uniform component L appears in the field direction. After integration over L the action llKit acquires 
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additional terms proportional to H, which can be recast into a shift of the time-derivative of n, as expected from the 
spins precession around the applied field: 



5(H) = S (d T n d T n + iH x n) = S + — [ [2iH • (n x d T n) - H 2 + (H • n) 2 l 

%9 C J 



(18) 



Observe that in the NLSM the constraint n 2 = 1 allows one to rewrite the last two terms of Eq. iflSl) also as H 2 nj_, 
where is the component of the order parameter perpendicular to the field. 

The non-linearity of the model l|16ll - ljl8Jl resides in the constraint n 2 = 1 for the staggered field. We will implement 
it by means of a Lagrange multiplier A(x, r), which corresponds to add a term J i\(n 2 — 1) to the action l|T8jl . and 
to perform an additional functional integration over A in the partition function^ In the saddle-point approximation 
A can be taken as a constant Ao, and one can expand the field n in terms of fluctuations around a given equilibrium 
configuration no. Both the value of the Lagrange multiplier and of the order parameter at each temperature will 
be determined by minimizing the action with respect to them. The result follows straighforwardly in the case of 
no magnetic field: assuming no = (0, cro,0), and integrating out in momentum space the (transverse) Gaussian 
fluctuations around it, n = (n a , ao, n c ), we obtain Z = exp(— <S) with: 

S= 1 -TrlogG+^ c [m 2 (al-l)], (19) 

where A is the area of the 2D system, m 2 = iX^lgc and the matrix G _1 is given by: 

I f ul+ c 2 k 2 + ml +m 2 



G - ~ I uj 2 + c 2 k 2 + m 2 + m 2 I ' (20) 



where w„ = 2tttiT and q are the Matsubara frequencies and the momenta, respectively, and the trace in Eq. l(T9ll is 
over w m ,q and the matrix indexes. By minimizing the action ijl9JI we obtain two equations: 



m 2 (To = 0, 

4 = 1-NI X , (21) 
where I± = (1/2)(J Q + I c ) accounts for the transverse fluctuations, with 

^ GaAK^n) = ^(K c (k,^)| 2 >, (22) 

where G a , c = (« 2 . c ) i s the Green function for the a, c mode. From Eq. l|2"Tll we see that two regimes are possiblei 22 i 2 ^ 
(i) uo = 0,™/0, which corresponds to the paramegnetic phase. Here m 2 plays the role of the inverse correlation 
length, defined by the second of Eqs. ll2Tll for 00 = 0, i.e. 1 = NI±(rn 2 ); (ii) m = 0, ao ^ which is the ordered phase, 
where the order parameter is an increasing function of temperature below T/v, defined as the temperature at which 
the mass first vanishes, i.e. 1 = NI±(0). Observe that in the 2D case the functions I a , c can be evaluated analytically 
and are given by: 



gT ^ f sinh(cA/2T) \ 
2?rc n \ sinh(w ajC /2T) J 



(23) 



As far as the magnon gaps are concerned, they are defined as the poles of the spectral function at zero momentum: 

Aa,c(w) = ImG a Jiu) n —> u + i0 + , k = 0) = [S(u - m a , c ) - S(uj + m a<c )] , (24) 

7r 2m a , c 

where the last equality only holds in the case of a matrix for the transverse fluctuations having the diagonal structure 
of Eq. il2f)ll . Thus, one can identify m a , c as the magnon gaps at zero field. 



A. Transverse field 



Let us analyze how the previous results are modified in the presence of a finite magnetic field. We first consider the 
case of a transerse field, for example H || a. Eq. Q18J1 then reads: 

S = So + — [ \2iH(n b d T n c - n c d T n b ) - H 2 + H 2 n 2 + i\2gc{n 2 - 1)1 . 
2ffcJ 
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As a consequence, using again n = (n , erg, n c), the only modification to the previous set of equations is the replacement 
of m 2 with m 2 + H 2 in the Eq. ll20l) defining the transverse fluctuations. Thus, one recovers the same phase transition 
as before (with negligible quantitative corrections to the Ijv and the cr value). As far as the magnon gaps are 
concerned, we see that the c mode is unchanged, while it occurs a shift of the mass of the a mode, leading to two 
poles at: 

2 2 i 7"T"2 2 2 

u a = m a + H , lu c = m c , 
corresponding to the result © that we derived above. 

B. Longitudinal field 

In the case instead of a longitudinal field, H || 6, Eq. l|T8jl reads: 

S = S + / [2iH{n a d T n c ~ n c d T n a ) - H 2 {n 2 a + n 2 c ) + iX2gc{n 2 - 1)1 . 

Thus, when n = (n a , <jq, n c ) we find that the equivalent of the inverse matrix 12011 acquires off-diagonal terms propor- 
tional to the applied field: 



G 



-l 



1 f c 2 k 2 + ml - H 2 +m 2 2wH 



gc 



-2loH uj 2 + c 2 k 2 +m 2 -H 



2 



As a consequence, the matrix G which defines the transverse fluctuations reads: 

where e 2 c (k) = c 2 k 2 + m 2 c . Due to this structure, the poles of the spectral functions for the a, c modes are the zeros 
of the determinant of the G matrix at k = and iu> n — ► uj + i0 + : 

{-to 2 + e 2 (0) - H 2 ){~lu 2 + e 2 (0) - H 2 ) - Alu 2 H 2 = 0, 

and correspond to the two solutions (1 1 Of) determined above using semiclassical spin-wave theory. Observe that in 
principle both solutions appear in the spectral function of the a or c mode. However, the spectral weight associated 
to the two solutions uj a ,c differs in the two cases. For example, for the a mode we have: 



Aa{iO > 0) 



-^-5(UJ - UJ a ) + 7T^-S(U ~ 0J C ) 

2u a 2uj c 



(26) 



where the residua at the poles are (see also Eq. 14111 below) Z a c = ±(— uj 2 c +m 2 —H 2 )/(u; 2 —oj 2 ). Since Z a /uj a Z c /w c 
one can conclude that the spectral function of the a mode is dominated by the pole at u a , and conversely for the c 
mode, confirming the identification of the two function iflOj) as the correct magnon gaps in an applied longitudinal 
field. 

C. In-plane field 

Finally, for the sake of completeness we analyze the case when the field is applied in the plane forming an arbitrary 
angle <fi with the a axis. Thus, the field has both a longitudinal (H sin 4>) and a transverse (H cos</>) component, and 
we expect an intermediate behavior between the two cases analyzed above. Following the same line of calculculations 
described in the previous subsections, we obtain: 

, m 2 , + m 2 H 2 , 



H 2 H 4 

—A(m\ +m 2 ) + m 2 L H 2 sin 2 (j) - m 2 H 2 cos 2(f) + — [A 2 + B] , (27) 
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where 

A = 3 sin 2 cf> + cos 2</>, 
B = A cos 2<p sin 2 0. 

For the c mode we obtain an analogous expression, with a plus sign in front of the square- root term in Eq. 1I27JI . Thus, 
Eq. ll2"Tj) reduces to the results © and ijlOj) when = and <j) — tt/2, respectively. The behavior of cu a for various 
values of the angle as a function of the field strength H is plotted in Fig. [21 At 4> — tt/2 the expression I|27J1 (and 
then also Eq. lUOJl l is vanishing at a field H c — m a . Indeed, as we discuss in detail in Sec. V-A, at this critical field 
the spins perform an in-plane spin-flop transition to orient perpendicularly to the magnetic field. When 4> deviates 
from 7r/2 the longitudinal field component decreases and the spin-flop transition moves to a higher value of the field. 
Accordingly, the field dependence of the gap changes continuosly, going from a decreasing function to an increasing 
one, recovering at = the increasing behavior dictated by Eq. 10, characteristic of a purely transverse field. From 
Fig. 12 it is also clear that the two extreme cases are also the ones where the largest deviation of the gap from the 
zero-field value can be observed. 

0.04 
0.035 



0.03 
0.025 

0.015 
0.01 
0.005 


5 10 15 20 25 30 

H(T) 

FIG. 2: (Color online) Field dependence of the in-plane magnon gap for an in-plane field applied at an arbitrary angle (f> with 
the a axis, see Eq. 1271 . The case <f> = corresponds to the transverse-field case analyzed in Sec. IV- A, Eq. 0, while the case 
4> = 7r/2 corresponds to a purely longitudinal field, Eq. which has a spin-flop transition at H — m a . The behavior of the 

gap above the spin-flop transition is described in Sec. V-A below. 




V. NLSM WITH INTERLAYER COUPLING AND DM INTERACTION 

Once we have established the equivalence between the semiclassical approach and the NLSM derivation of the spin- 
wave gaps for an ordinary easy- axis antiferromagnet, we can discuss the more general case of the model © where 
also interlaye r coupling and DM interactions are present. The NLSM for this system has been already derived in 
Refs 1 1(1 121231 where it has been shown that in the absence of magnetic field a 3D analogous of Eq. Illlit holds: 



So = ±- J2 [ {(dr^mf + c 2 (Vn m ) 2 + r?(n m - n„ l+1 ) 2 + {D + n a m f + T C «J 2 } , 



(28) 



where m in the layer index, D + = £> + x a , D + = ASd + , T c = 32JS 2 (Ti — T^) and r\ = 2JJ±. In this notation the 
in-plane and out-of-plane modes have masses: 



m a = D+ = ASd +1 m c = y/T c = ASy/2J(T 1 - T 3 ), (29) 



in agreement with the results lTL2]) above. Moreover, the uniform magnetization of each layer acquires an additional 
contribution proportional to the DM vector D + with respect to Eq. I|17J1 : 

F (n m xri m ) + -^[(-l) m n m xD++H-n m (n m -H)], (30) 



8aSJ "* 8aSJ" 
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where the oscillating factor (— l)™ 1 in the term proportional to D + accounts for the effect of the tilting of the ochtaedra 
on neighboring planes, as discussed below Eq. <EJ|. As one can easily see applying the saddle-point approach described 
in the previous Section to the action J2HJ, at H = the system orders AF below T/v in a 3D staggered configuration, 
with n m along b in each layer. Moreover, due to the oscillating factor (— l) m in Eq. l-'>()ll . the spins in each layer 
acquire a FM components L m , with the vectors L m ordered AF in neighboring layers, see Fig. 03 

H-0, <„„> = (L m , = (-i r (^± = (-ir£g±. ( 31 , 

The additional term in D + in Eq. ll-'iOt translates in an additional coupling between the order parameter and H when 
the full NLSM action at finite magnetic field is computed: 



5(B) = So + ^- / [ 2lU ■ ( n ™ x 9 - n '") - H 2 + (H • n m ) 2 - (-l) m 2H • (n m x D+)] . 



(32) 



As it has been discussed in Refs. Isfll0ll2ll4 the effective stag gered field acting on the AF order parameter due to 
the presence of the DM interaction makes the system an unconventional easy-axis AF. As far as the spin- waves gaps 
are concerned, the results of the previous Sections apply only in some specific regimes, as we shall analyze below. It 
is worth noting that the last three terms in Eq. 1(32)1 are proportional to: 

-H.(l£+L™), (33) 

where and L^ M are the contributions to L m in Eq. ll-'iOt proportional to the magnetic field and to the DM term, 
respectively. As a consequence, the ground-state of the action l(32|l will be determined by the competition between 
the energetic cost of the bare action So and the tendency of the system to maximize the uniform spin components in 
the field direction, to gai n energy f rom the term ll.'i-'ill . Even though part of the ground-state phenomenology has been 
already described in Ref. l8lflUlCMiaTl4L here we will derive these results within the general language of the saddle-point 
approximation for the NLSM, by focusing on the magnon-gaps behavior^ that will be compared with the expectation 
for an ordinary easy-axis AF, described in the previous Section. Then we shall also compute the effect of quantum 
and thermal corrections, which allows us to investigate the (H,Tn) phase diagram in the various case. 



b 

— > 



FIG. 3: (Color online) Spin configuration at zero applied magnetic field in La2Cu04- The arrows with a solid tip represent 
the staggered spin components ( — l) m e lCJ ' ri n™, the arrows with a two-lines tip represent the uniform spin components L™, the 
arrows with an open tip the full spins S™. 



A. H paralell to b 

When the field is in the b direction the last term in Eq. (l32l) generates a staggered field along the c direction which 
leads to a rotation of the order parameter in the be planei&Ifiiii 

5(B) = So + ±- ]T f [2iH{n a m d T n c m - n c m d T n a J - H 2 [(n c m ) 2 + «) 2 ] + 2h m rf m + a m 2 3 c(n 2 „ - 1)] , 

where we put h m = (~l) m HD + and we introduced explicitly also a set of Lagrange multipliers X m which enforce the 
constraint n 2 ^ = 1 in each layer. Due to the anomalous coupling, a finite component in the c direction can arise. For 
a generic configuration n m = (0, a h m , <t£J we find the ground-state action (m 2 ^ = 2i\ m gc): 

= ^ E - ff ™+i) 2 + ^< " CT »+i) 2 ~ 2hm < + (™ 2 ~ # 2 )«J 2 + mUM 2 + «J 2 - 1]. 
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The saddle-point equations then read: 

?7( 2cr m ~ ~ CT m-l) + = 0, 

(ml - /f 2 )< rl + /i m + V (2a^ - <„ +1 - c^.J + m 2 <„ = 0, 

(^) 2 + KJ 2 = i- 

At low field, one can easily check that the classical configuration is given by an order parameter with a uniform o\> 
component in neighboring layers and an oscillating a c component (which corresponds to the c components of the spins 
coming from the (— l) m e lC *' ri n™ term ordered ferromagnetically in neighboring planes, see Fig. 

Sd = -7T-W% m 2 n = m 2 , 
2gc 

tt ry 

= = l-o*, a c m = (-l) m a c = (-ir + (34) 

In this configuration the L m vectors are given by: 

L ° M = ^h [K]5lb + hir(Tb5lcl (35) 

so that the average (i.e. summed over neighboring layers) uniform magnetization is along b and given by (L) = 
(l/8aSJ)[D + \a c \ + 7fcr 2 ]x{,. Thus, the oscillating character of tr^ allows for the DM-induced magnetization to allign 
in the direction of the field, see Fig.|3 After inclusion of the transverse a, c fluctuations, the order-parameter equations 
read: 



a 2 = l-a 2 -NI A _{m = 0), 

T<T N , (36) 



HD+ 



m 2 + 4i] - H 2 ' 
below TV and 

<r b = (),=>• 1 = al+NIj_(m), 

T>T Nl (37) 



HD+ 



ml + 4?y — H 2 + m 2 ' 

above TV, where I± is computed using the matrix (12ofl discussed above for the case of longitudinal field. Moreover, 
since the system is now 3D, we have that the energy dispersion of the transverse modes is 

e 2 c (k, kx) = c 2 k 2 + 2ry(l - cos kj_d) + m 2 acl (38) 

where d is the interlayer sapcing and an additional integration over out-of-plane momentum k± must be included in 
computing I a<c . Thus, taking into account the non-diagonal character of the fluctuations matrix QSJ, we have for 
example for the a mode: 

T 1 /, / i , m2\ 1 Z a (k,k±) (3uj a (k,kj_) Z c (k,k ± ) j3w c (k,k±) 

Ia = ^T7 1^ Kk.Mi = tt V o — 7z u \ coth o + o tcnr~\ coth 5 > 39 

0V u £ k± y k ^2^ a (k,fc ± ) 2 2u, c (k,fc ± ) 2 

where V is the 3D volume of the system. Here w 0;C are the generalization of Eq. IIOII at finite momentum: 



£ 2 +£ 2 //c2_J\» 



^ iC (k,fcx) = ^ T ^+^±y^^ T ^j +AH 2 \^^y (40) 

where the explicit dependence of e a , c on momenta has been omitted. Analogously, the spectral weights Z a , c of the 
two poles are given by: 

-oj 2 ac + el(k,k ± ) -H 2 

Z , c (k, k±) = ± r „. r , (41) 

w 2 (k,fc ± ) - w 2 (k,fc ± ) 
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FIG. 4: (Color online) Momentum dependence of the spectral weight corresponding to the poles at w a (k,fcx) or u c (k, k±) in 
Eq. 13911 . Here k± = 0, m a = 0.02 J, m c = 0.05J, c = J, and we used units such that H=l T corresponds to 10 _3 J. Inset: 
momentum dependence of Z a ,Z c defined in Eq. 14111 . Observe that even thought Z a ~ Z c at H = 10 T as ka ~ 0.5 the 
contribution at k = is always predominant in the momentum sum 1391 . due to the fact that the a gap softens as the field 
strenght increases. 



and are plotted in the inset of Fig. 21 for k± — as a function of ka. Note that since in all the formulas the magnetic 
field is measured in units of g s ^B ~ 0.1 meV, and since typical values of J are of the order of 130 meV, we used for 
simplicity the equivalence H — 1 T=10~ 3 J. Due to the thermal factors, the main contribution in the momentum 
integration in Eq. I|39|) comes from k = k± = 0, where also the factors Z a /oj a have the largest value, see Fig. 31 As a 
consequence, we can safely approximate the momentum dependence of Eq. Il40fl with 

u£ c (k, k±) « u 2 a c + c 2 k 2 + 277(1 - cos k±d), (42) 

where uu a ^ c are the magnon gaps given in Eq. ijlOjl . and we can neglect the momentum dependence of Z a ^ c in Eq. fjj2j. 
We thus obtain below T/y: 

la = Z a 

(0)J 3 £>(Wa) + Z c (0)I 3D {u c ), (43) 
where is the extension to the layered 3D case of the integral (J2 



hn{M) = f- r % ml - Sinh(cA/2T) . (44) 

2ttc 7_ 7r 2tt I shih ^M 2 + 2tj(1 - cos z)/2T f 

Above Tjv we simply substitute lo\ c —> uj\ c + m 2 in Eq. 14-'tt . where as usual m 2 plays the role of the inverse 
correlation length, to be obtained solving the self-consistency equation Q37JI . As far as the c fluctuations are concerned, 
the previous result is clearly reversed, with a larger contribution coming from the pole at oj c , since now Z a c (0) = 
T(-< c + m 2 -ff 2 )/(a, c 2 -c 2 ). 

As the field strength increases we see that, according to Eq. IIOII . the smaller gap ui a decreases, and vanishes at the 
critical field: 

HP=m a = D+. (45) 

Indeed, at we have a spin-flop transition: the in-plane component of the order parameter rotate from the b to 
the a direction, so that the classical configuration becomes: 

n m = (a o ,0,(-l) m cr c ). (46) 

The uniform magnetization changes correspondingly: 



13 



0.005 



0.004 



0.003 

A 
_l 

V 

0.002 



0.001 





5 10 15 20 25 30 35 

H (T) 

FIG. 5: (Color online) Magnetic-field dependence at of the spins ground state and of the T — average uniform magnetization 
in the field direction for a field along the b axis. At finite magnetic field the order-parameter n rotates in the be plane. The 
spin configuration is determined as usual by Si/S = ( — l) m e lQ ' ri n™ + L™, where the staggered and uniform components are 
indicated by arrows with different tips as in Fig. the arrows with a two-lines tip represent the uniform components L m and 
the arrows with a solid tip the staggered components ( — l) m e lC! ' ri n* n . At H < » 18 T (left in the figure) n m is defined in 
Eq. 134t and L m in Eq. 1351 . so that the spins have both the staggered and uniform component in the be plane. The uniform 
components are due both to the magnetic field and to the DM interaction, see Eq. 13511 . and the sum of Lai +L m+ i is alligned 
along the field. Above the spin-flop transition (right in the figure) the staggered components n m lie in the ac plane see Eq. 
1461 (so that the part along a is orthogonal to the plane of the figure and indicated by a cross) and the uniform components 
L m point along 6, see Eq. H47fl . The order-parameter values used to evaluate the (L) are the same reported in Fig.Q Observe 
that since in our approximation the critical field He is independent of temperature, as well as cr c , the uniform magnetization 
is temperature independent below Tat. 



Uniform magnetization 




and jumps discontinously at the spin-flop transition by a quantity H[l — o c {H^)]/&aSJ ', see Fig. El 

By adding b and c fluctuations around this ground-state solution we obtain the new saddle-point equations: 

a 2 a = l-a%- Nlj_{m = 0), 

a, T<T N , (48) 



and 



4r? 



<J a = 0,=M = a 2 c +NI x (m), 

°c = 2^ A HD+ 2^ 2 ' T>T N , (49) 

where I± accounts for the b, c fluctuations described by the inverse matrix: 

A_i _ ( u 2 {q)+H 2 -ml +m 2 \ , , 

\ w 2 (q) + m 2 c -m 2 a +m 2 )> [b[j) 

where w 2 (q) = w 2 +c 2 k 2 + 2r;(l — cosk±d). As far as the magnon gaps are concerned, we see that above the transition 
the field in the b direction acts as a transverse field, since the direction of the magnetization has changed. Once again, 
the behavior of the magnon gaps can be easily read from the Green's function matrix We find that the in-plane 
mode corresponds now to a fluctuation of the 6 component, with a field-dependent mass, while the out-of-plane mode 
does not depend on the field but should be rescaled with respect to the m a gap: 



= uj 2 = y/lP-ml, lu 2 = y/m* - m 2 . (51) 



The resulting field dependence of the magnon gaps is reported in Fig. El 
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FIG. 6: (Color online) Field dependence of the magnon gaps at T = for a field applied along the b axis. We used m a = 
0.02J, m c = 0.05J, and units such that H=l T corresponds to 10 J, so that the first critical field He = D+ — m a = 20 T. 
Below He the in-plane and out-of-plane mode are correspond to a and c component fluctuations, whose gaps are given by ui a 
and uj c in Eq. 1101 . respectively. Above the critical field the in-plane mode corresponds to fluctuations of the b component, and 
the magnon gaps are defined by Eq. 15111 . 

By further increasing H one will finally reaches a second critical field above which the transverse in-plane component 
of the staggered order parameter vanishes, a a — 0, but the <r c component is still nonzero, see Fig. Since by definition 
|n| < 1, we can see from Eq. (14811 that the solution Il41it is valid only when a c < 1, so at fields lower than&ifi: 

H?>™™) = {ml + A V -ml)/D + . (52) 

However, the estimate (§2) does not take into account quantum fluctuations, which reduce the T — value of the 
in-plane order parameter according to Eq. !■' >(ill . Since the transverse fluctuations I± do not depend strongly on the 
magnetic field, one can see that the second critical field, defined as the field at which a a (T — 0) = in Eq. H48JI . is 
given approximately by: 

^ 2 ^ mg+ ^' m V ( g ^O,T^O), (53) 



where <Jo(H = 0) = \/l — N I± (T — 0, H = 0) is the magnetization of the system measured at T — along the b 
axis without external magnetic field. Observe that quantum corrections can indeed reduce considerably the second 
critical field. In Fig. we report as an example the phase diagram for La2Cu04 of the transition temperature T/v vs 
H obtained in the saddle-point approximation. By estimating the parameter values from the Raman measurements 
of the magnon gaps in Ref. [l3 (see alsoGU, we choose m a = 0.0116 J, m c = 0.034 J, rj — 8 x 10 5 J 2 , with J — 130 
meV and g s = 2.1, as appropriate for the b direction. For the stiffness and spin-wave velocity we use p s — 0.07J and 
c = 1.3J respectively, which are not far from the standard values quoted in the literatur o 2 ^ 2 ^ and have been shown 
to be appropriate to reproduce (in the same approximation) the uniform-susceptibility datai^. In the inset we also 

report the T = value of the order-parameter components as a function of the field. As we can see, below H^ l) =D + , 
a a = and er^ ^ 0, while the situation is reversed above the spin-flop transition. The component a c is continuous at 
the spin-flop transition. In principle, its slope as a function of the magnetic field changes at the spin-flop transition 
according to Eqs. ll36|l and l|4*8)l. However, for the parameter values used here, as appropriate for La2Cu04, this 
change is almost undistinguishable in Fig. [3 Moreover, also the magnitude of the in-plane component is continuous 
at the transition, the change being only in its direction. 

As far as the second critical field is concerned, it turns out that Eq. (1531) . which uses the value of oq at H — 0, is an 
excellent estimate of the second critical field He , obtained by means of the self-consistent value 1 — NI±(H): indeed, 
since we found cr (T = 0, H = 0) = 0.4 (see inset), and H { ?' nalve) = 77 T, Eq. ® would give H^ 2) = 30.8 T, which is 
almost the value found numerically, see Fig. [3 It is worth noting that has been measured recently in 1% doped 
La2- a; Sr 3 ;Cu04 sample by Ono et ali2&, who found — 20 T. Thus, even though has not been measured in 
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FIG. 7: (Color online) Phase diagram of the system for H || b. The parameter values are chosen to reproduce the data of 
La2Cu04, i.e. D+ = m a = 0.0116J, m c = 0.034J, J = 130 meV, p s = 0.07J and c = 1.3Ji&. For the interlayer coupling we 
used r\ — 8 x 10~ 5 J 2 , which correctly reproduce the spin-flop transition measured in Ref. 15 for H || c, see end of Sec. V-B. The 
magnetic field is converted in energy trough H — > g s ^BH, where fiB is the Bohr magneton and g 3 = 2.1 is the gyromagnetic 
ratio in the 6 direction. Inset: field evolution of the staggered order parameter components at T — in the various phases. 



undoped La2Cu04, where it is expected to be larger than in the 1% doped sample, the value of 30 T following from 
Eq. H53H seems more realistic than the bare estimate 15211 . which gives 77 T. Observe that neglecting the quantum 
renormalization of the order parameter in the estimate of the second critical field can lead to an underestimate of the 
mass of the c mode, as it has been done in Ref. |J where Eq. l(52|l has been used. 

Fibally, we note that in the saddle-point approximation used so far the transverse gaps are constant in temperature 
below Tpf. However, one would expect that a better approximation could reproduce the softening of the transverse 
gaps as the temperature increases, as observed experimentally. In this case also the value of the first critical field 
He 1 - 1 would acquire a temperature dependence, which is instead absent in the phase diagram of Fig. Moreover, this 
could also smoothen the discontinuity of Tjv at the spin-flop transition found at saddle-point level. 



B. H parallel to c 



When H is along the x c direction the last term in the action Ij32l) gives rise to an effective longitudinal staggered 
field: 

5(B) = So + -L I [^H(n b m d T n a m - n a m d T n h J - H 2 + i/ 2 «J 2 ] + 2h m n b m + *A m 2 5 c(n 2 rt - 1)] . (54) 

Following the line of the analysis performed in the previous section, we first determine the ground-state configuration 
in the presence of the magnetic field. Since also the effective staggered field is longitudinal, we do not expect in this 
case to have a change of direction of the equilibrium configuration. We can then look for a ground-state solution of 
the form n r „ = of n X6, which gives: 

S * = |^ E ^ ^+i) 2 + 2^ m < + mU(v° m ) 2 - 1]. (55) 

The saddle-point equations then read 

r,(2a° m - a° m+1 - ^_ 1 ) + h m + = 0, 

«) 2 = 1. (56) 

As a consequence, two solutions are possible: (i) the order parameter is the same in all the layers, and the ground-state 
action S c i vanishes: 

al = a = -l, m 2 = -(-l) m h/a = {-l) m h, S cl = 0. (57) 
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This configuration is the same of the case without magnetic field. The uniform spin components are: 



H 



H 



DM 



8aSJ 



(-iy 



8aSJ 



Fo|x c 



(58) 



so that the L^ M , which are ordered antiferromagnetically in neighboring layers (see Fig. 03), do not contribute to the 
average uniform magnetization (L) = h H (see Fig.JHJ) , but one takes advantage from the out-of-plane antiferromagnetic 
coupling; (ii) the order parameter change sign in neighboring layers, which means that the spins order ferromagnetically 
in the c direction, with the moment L^J M oriented in the same direction: 



H 



H 

8aSf 



DM 



Da 



8aSJ 



|oo|x c 



(59) 



giving (L) = {H + D+\cr \)/8aSJ. This spin flop of the uniform c components of the spins leads to a lowering of the 
energy when the gain in magnetic energy is larger than the cost coming from the interlayer AF coupling. Observe 
that the average uniform magnetization jumps discontinuosly at the spin-flop transition, the jump being proportional 
to (To, so that it decreases as the temperature increases, see Fig. 03 The classical action in this configuration is: 
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FIG. 8: (Color online) Magnetic-field dependence of the spin configuration and of the average uniform magnetization at two 
different temperatures for a field along the c axis. The spin configuration is determined as usual by Si/S = (— l) m e lCJ ' ri n m +L m , 
where the staggered and uniform components are indicated by arrows with different tips as in Fig. and Fig. At H < H c , 
where H c = 6 T at T = 10 K and H c — 4.6 at T = 200 K, see also Fig. llOl below. the system has 3D antiferromagnetic order as 
described by Eq. 15711 . The difference between the L m components (arrows with a two-lines tip) in neighboring planes is due 
to the DM-induced term L m M which changes sign from one layer to the next one, while is always parallel to H, see Eq. 
JHHJ. Ab ove the spin flop the spins are ordered ferromagnetically in neighboring layers, see Eq. 1GIH . allowing for the uniform 
L m components to allign along the field in all the layers. The jump at the transition is proportional to the order parameter 
o"o, and decreases as the temperature increases, see inset of Fig. 1101 



a° m = (-l) m a = -1, m 2 m = m 2 = -h/a - 4 V = h - Ar,, S cl = (-h + 2r?) , (60) 

where iVj is the number of layers. When h > 2r\ this second solution becomes energetically favorable, so that the 
critical field for this spin-flop transition is: 

Hc = ^-- (6i) 

When transverse fluctuations are included the first of Eqs. 15 (ill will not change, while the value of the order parameter 
do will acquire quantum and thermal corrections due to transverse fluctuations, according to: 

a 2 = l-I ± {h/<j ). (62) 
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In Eq. I(i2t we included the explicit dependence of the transverse fluctuations on the order parameter oo via the 
effective staggered field h, as we will derive in detail below. To determine I± we should distinguish the case below 
and above the spin-flop transition. Let us first consider the case ifS?)) of low field. Since the Lagrange multiplier is 
oscillating in neighboring layers, it gives rise to a coupling between the transverse modes at k± and k± + Q±, where 
Q± = ir/d. We then have: 

i 

S = Sd + - £ ^*i(k,fe x )G- 1 * a (k,fcx), 

a— a,c k,k± 

where ^(k, k±) = ((n a (k, k±), n a {\s., k±_ + Q±)), the sum J^k 1S over the reduced Brilluoin zone for k±, i.e < 
k± < ir/d. The Green's function for each mode is now a matrix. For the a mode we find: 

1 / c 2 +c 2 (k,fcj 



with w 2 (k, fc_i_) = to 2 + c 2 k 2 + 2r;(l — cosk±d) and m 2 = —h/ao, according to Eq. {SJJ. The inverse Green's function 
for the c mode has an analogous expression, except that m 2 — > m 2 + H 2 . After integration of the Gaussian fluctuations 
one obtains the order-parameter equation 162L provided that Eq. iKi-'it is used to compute the fluctuations: 

(|n a (w„,k,fc ± )| ) = — j- . (64) 

det G a 

The integral of the a- mode fluctuations has a structure similar to the one of Eq. ftfflU • with two contributions at the 
eigenvalues of the matrix G a ■ 

1 \ -> „ , , , .,9, 1 \ Z + (k±) , Bu)+(k, k±) Z-(k±) , /3w_(k, . , 

J " = ^ E d««(^ k .yi 2 ) = ^ E 4^) coth 2 + Mi^ c ° th ^' (65) 



where 



^2 (kfc±) = £g(k,fc ± )+£ 2 (k,^+Q ± ) ± // £ g(k,fc ± )-eg(k,fc ± +Q ± ) 



= c 2 k 2 + m 2 + 2?7 ± cos 2 k ± d + m 2 . (66) 

Analogously to the case discussed in the previous section, the spectral weights Z± of the two poles are not equivalent, 
and determine the main character of the excitation. In this case we have: 



, , _ -QJ± +^ 2 (k, k± + Q±_) _ 2?/ cos k^d =p \J W' cos 2 k±d + m 2 
£ ^ ~ T w 2 (k,/c ± )-w 2 (k,fc ± ) ~ T cos 2 + to 2 ' 

The momentum dependence of Z± is reported in the top panel of Fig. EJ while in the lower panel the two solutions 
uj±(k = 0, k±) are plotted. As one can see, as k± increases the spectral weigh of the momentum sum in Eq. j^H 
moves from the solution u>-(k±) to the solution ui+(k±), which follow closely the bare function oj(k±) in the two 
regimes < k±d < ir/2 and ir/2 < k±d < ir respectively. The effect of the magnetic field is then twofold: it affects 
the magnon gap k± = and opens an additional one at k±d = ir/2. To compute explicitly the momentum sum in Eq. 
(16511 we can observe that Z±(kj_) only depend on the out-of-plane momentum k±. We can then perform the usual 
integration over the in-plane momentum k in Eq. I|65J1 . obtaining an expression similar to Eq. 143|l - 144t : 



where 



gT r dz „ , . , f sinh(cA/2T) 1 _ . . , f sinh(cA/2T) 1 TT TT 



A± (z) = yjml + 2r] ± y/4r] 2 cos z + (HD+/a ) 2 . (69) 
For the c fluctuations one finds the same results, provided that m 2 — > m 2 + if 2 in all the above formulas. 
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FIG. 9: (Color online) Low panel: out-of-plane momentum k z = k± dependence of the two solutions w± of Eq. Utitill at k = 0, 
compared to the momentum dependence of the function u>(k = 0, k±) = \J m% + 27/(1 — cos k±d). We used D+ = m a = 0.02J 
and H — 1 T= 10~ 3 J. Observe the small difference of lu- with respect to uj at k± = 0. Top panel: momentum dependence 
of the spectral weight Z±(k±) defined in Eq. H(i7ll . As one can see, at small momentum Z- 3> Z+, while at k±d > ir/2 the 
situation is reversed. As a consequence, the momentum sum in Eq. Ifi'it select always between u)+ and w_ the solution which 
follows most closely the momentum dependence of the solution uj(k±) of the case H = 0. 



Let us discuss now the issue of the magnon gaps. The spectral function of the a mode at k = 0, fcj_ = has a 
two-poles structure analogous to Eq. l|2fijl . i.e.: 

a , b n \Z+{k^=0) sf , Z_(fc ± =0) r . 
A ac (uj>0)= y j -S(w-oj + )-\ K — '-8{w-u-) 

However, as observed before and shown in Fig. at fcj^ = is Z_/cj_ ^> Z + /lo + , so that only the second term 
contributes in the previous equation and allows us to identify the magnon gap as the k = 0, k± = limit of the the 
function w-(k, k±) above, i.e.: 



2r) — \ 4r] 2 



H 2 + 277- W Arf 



H<H C , 



(70) 



which reduce to the conventional ones when H — 0. Here we used explicitly that m 2 = —h/ao, according to Eq. {29 < 
Observe that this result is quite different from the interpretation given in Ref. Eli, where it was claimed that the 
acoustic and optical modes are mixed. Instead, at k± = only the mode w_ is observed, as Raman measurements 
confirm Moreover, we stress that according to the discussion below Eq. 10, for an ordinary easy-axis AF we 
expected that uj a is unchanged and ui c hardens for a field parallel to c. Instead, due to the presence of the DM 
interaction, the two modes have a field dependence ~ — -ff 2 7a,e, with j a — D^/^a 2 ) and -f c = — 1 + / (4t](Tq) . As 
a consequence, the a mode always softens, while the behavior of the c mode depends on the ratio D+/ '(Ai]a ). 

Let us analyze now the case H > H c . According to Eq. KiOII we find a uniform saddle-point value of the constraint, 
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T(K) 



FIG. 10: (Color online) H c — T phase diagram for H || c. We used the same parameter values as Fig.0 except that in this case 
g s = 2.4, as appropriate for a field in the c direction-^. The order parameter ao(T) has been determined self-consistently trough 
Eq. 1621 . computing the transverse fluctuations with Eq. I(i8ll and 173fl below and above the spin-flop transition, respectively, 
where the critical field H c is defined trough Eq. 17411 . Inset: temperature dependence of the staggered order parameter at 
ff = lT. 



?7i„ t = m 2 = —h/ao — 4r], so that the matrix of the transverse fluctuations is again diagonal in momentum space: 

r-i _ { ul + c 2 k 2 -2i 1 (l + cosk ± d) + m 2 a -h/cro \ , > 

\ col + c 2 k 2 - 2ry(l + cos k±d) + m\ + H 2 - h/a )' [ ' 

However, in the classical configuration the spins are order ferromagnetically in neighboring planes, see Fig. IH1 this 
means that the low-energy spin fluctuations are those at k± — ir/d in the notation of Eq. itTTll . so that the spin- wave 
gaps evolve at H > H c in: 



2 2 

uj„ = m„ 



( HP, 



u 2 c = m 2 c +H 2 +(^±y H>H C . (72) 

Since the matrix H71H admits two simple poles, the transverse fluctuations are described by the function 14411 . whith 
the masses given by the previous equation: 

I a ,c = hD(u a ,c), H>H C (73) 

Observe that both for the case H < H c and H > H c the integral of transverse fluctuations, given by Eqs. l(68jl and 
173|l respectively, depend explicitly on the Lagrange multiplier h/ao, so that Eq. I(i2ll is a self-consistency equation for 
the order parameter at all the temperatures. This is quite different with respect to all the cases analyzed before, where 
J_l depends only on the transverse masses and two separate regimes exist, m 2 = 0, ao ^ below T/v, and m 2 ^ and 
(To = above Tjv (see Eq. ^211 ). Here instead, due to the effective longitudinal field h m in Eq. (1o4t . instead of the two 
regimes we obtain a single self-consistent equation Illi2li valid at all the temperatures. As a consequence, exactly as for 
a ferromagnet in the presence of the external field, the order parameter never vanishes, and the transition transforms 
into a crossover from a regime where ao is large to one where a is smalU^. This is shown in the inset of Fig. El where 
we obtained ao(T) at H = 1 T by solving numerically the self-consistency equation 116211 at various temperature. It is 
worth noting that once that transverse fluctuations are included, also the definition of the critical field llfilll changes. 
Indeed, since in general ao is lower than 1 (including T = due to quantum correction), the critical field becomes 
itself a function of temperature. A first estimate of this effect can be done by evaluating again the value of the action 
S c i in Eq. KiOII for a generic ao- We then obtain that the critical field is: 

4?7 ap 

H -W + JiTaJr ^ 
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A more precise evaluation of H c could be done including also the contribution of the Gaussian transverse fluctuations 
to the action l(fi"3jl . However, already Eq. itTHl allows one to recognize that as the temperature increases the decrease 
of the order parameter ao induces also a decrease of the critical field for the spin-flop transition. The resulting H c — T 
phase diagram, obtained by means of Eq. I|74J1 where ao is the solution of the self-consistent Eq. Il(i2t . is reported in 
Fig. EH 

Once that we determined self-consistently the values of the order parameter and of the critical field, we can also 
compute the field dependence of the magnon gaps. In Fig. ifTTj) we show the field dependence at T = of ui a and u c , 
as given by Eqs. I|70|) and l|72|) below and above the spin-flop transition, respectively. For the interlayer coupling i] 
we used the value r\ = 8 x 10~ 5 J 2 , which allows us to obtain a critical field at low temperature around 6.5 T, as the 
one measured experimentally^. Observe that this value of r\ is quite similar to the one obtained in Ref. from the 
jump of the experimental measured in-plane gap at H c , even though such an estimate is done neglecting quantum 
corrections to the order parameter. With this parameter values we find that the below the spin-flop transition also 
the c mode is slightly decreasing. 
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PIG. 11: (Color online) Field dependence of the magnon gaps at T = for H \\ c. Here H c « 6.5T, so below it the gap decrease 
according to Eq. 17011 . while above it they increase according to Eq. I|72jl . The parameter values are the same of Fig. CSJ. 
Inset: field dependence of the zero-temperature order parameter. 



C. H parallel to a 

Finally, let us consider the case of a field along x a , i.e. parallel to the direction of the DM vector D+. In this case, 
the last term of Eq. l|32l) vanishes, and the system behaves as a conventional easy-axis AF. As a consequence, Eq. JjjJ 
holds, giving a hardening of the a gap and leaving the c gap unchanged. 

VI. CONCLUSIONS 

We have investigated the field dependence of the magnetic spectrum in anisotropic two-dimensional and 
Dzyaloshinskii-Moriya layered antiferromagnets. Starting from the appropriate spin-Hamiltonian for each case, we 
obtained the magnon gaps and spectral intensities as a function of the applied magnetic field, and we discussed the 
various possible ground state configurations and phase diagram. In particular, we showed that the peculiar coupling 
of the magnetic field with the staggered order parameter induced by the DM interaction gives rise to very interesting 
magnetic phenomena, such as spin-flop transitions and rotation of spin quantization basis. The predictions of the 
theory developed in this article are now ready to be compared with Raman spectroscopy experiments in Sr 2 Cu02Cl2 
and La2Cu04, reported in the forthcoming articled 
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